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ABSTRACT 


A frequency-hopping binary frequency shift keying (BSFK) ratio-statistic receiver 
with multihops per data bit is an effective electronic counter-countermeasures (ECCM) 
system against partial-band jamming interference. Interference 1s modeled as Gaussian 
noise. Orthogonal binary signaling and independent fading diversity is considered over 
frequency-nonselective, slow fading Rayleigh, Rician, and Gaussian channels. A forward 
error correcting coding scheme is implemented for a 1/2 rate convolutional code 
algorithm. The probability of bit error is examined for different levels of diversity, 
thermal noise, severity of fading, fraction of bandwidth jammed, and jamming power. 


Uncoded and coded system comparisons are done to determine worst case performance. 
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I. INTRODUCTION 


A frequency-hopping (FH) binary frequency-shift keying 
(BFSK) with multihops per bit is an effective electronic 
counter-countermeasures (ECCM) system. It provides low 
probability of intercept (LPI) and antijamming capability (AJ) 
for the communication channel. An effective strategy called 
partial-band noise jamming is used against this communication 
system [{Ref.1:pp.471]. 

The principal requirement of the communication system is 
to transmit the source information of a binary data sequence 
by means of M-ary FSK over the channel. We consider that the 
channel is modeled as a fading dispersive channel due to 
multipath problems and the jammer as an additive Gaussian 
noise. The addition of a frequency-hopping scheme provides 
some immunity against interference over an allowed system 
bandwidth of W Hz. 

Many different receiver structures for multihops per data 
bit of FH/MFSK waveforms transmitted over the channel have 
been considered in order to counter or mitigate the effects of 
partial-band jamming and fading dispersive channels. Examples 
of such receivers are: 1) square-law non-linear combining with 
adaptative gain control, 2) soft-limiter amplitude control, 


3) self-normalizing soft decision, and 4) ratio-statistic 


combiner. 

Keller and Pursley analyzed a ratio-statistic combining 
receiver ina frequency hopping spread spectrum system in the 
presence of partial-band interference for a channel with 
Rayleigh fading [Ref.2:p.145]. Their paper shows that the 
spread spectrum scheme provides some immunity to jamming and 
interference, and the ratio-statistic combining technique 
improved immunity to partial-band jamming. 

Riley analyzed the performance of a fast frequency hopped 
noncoherent BFSK receiver with ratio-statistic combining over 
a Rician fading channel with partial-band interference 
[Ref.3]. This thesis involves a comparison of the performance 
of envelope netectone versus square-law detectors, and the 
following conclusions were established: 1) ratio-statistic 
combining used in conjunction with diversity limits degrada- 
tion due to fading and partial-band interference and generally 
provides an improvement in overall performance, 2) Rician 
fading and partial-band interference have a significant impact 
on recelver performance, 3) ratio-statistic combining provides 
protection against partial-band interference as diversity is 
increased, and 4) for diversity of three and four, the 
receiver using envelope detection performs better than that of 
a receiver using square-law detection. 

Prior work by the authors [Ref.2}-[Ref.3] analyzed the 
performance of uncoded systems using soft decision receivers 


in the presence of fading channels with thermal noise and 


partial-band noise jamming. 

The purpose of this thesis is to show the performance of 
a forward error control (FEC) coded system uSing a ratio- 
statistic combining BFSK receiver with envelope or square-law 
detection in a fading channel with partial-band interference. 
The present work examines the probability of coded bit error 
versus different environment cases such as diversity, thermal 
noise, jamming duty factor, signal to jamming ratio, direct 
Signal to noise ratio, diffuse signal to noise, envelope and 
square-law detection. Performance curves are obtained for 


coded communication systems. 


II. COMMUNICATION SYSTEM DESCRIPTION 


A. MODEL. 

The general structure of the communication system to be 
evaluated is illustrated in Fig. 1.1. The system is charac- 
terized by an error probability which depends on the charac- 
teristics of the transmitter, the channel, and the receiver. 

This model is comprised of the following elements: 

e Sequence of binary messages. 
e Encoder. 

e M-ary FSK modulator. 

e Spread spectrum hopper. 

e Fading channel. 


e Partial-band interference modeled as additive Gaussian 
noise. 


e Thermal noise. 
e Spread spectrum dehopper. 


e M-ary FSK demodulator with square-law or envelope 
detection. 


e Ratio-statistic soft decision. 
®e® Decoder. 


e Error corrected binary sequence message. 


_ _ 


PARTIAL-BAND NOISE JAMMING —_ > Cc 
DECODED MESSAGE THERMAL NOISE a 





Figure 1.1 Communication System Model. 


B. PARTIAL“BAND INTERFERENCE. 

In an interference environment, frequency hopping spread 
Spectrum (FHSS) systems are used to reduce receiver perfor- 
mance degradation due to different type of jammers. In order 
to explain the partial-band interference term, it is necessary 
to understand the barrage noise jammer. This jammer transmits 
bandlimited Gaussian noise with a one sided power spectral 
density (PSD) of N; W/Hz and the jammer power spectrum covers 
exactly the frequency range of the spread spectrum signal. 
The effect of the barrage noise jammer on the system is to 
increase the Gaussian noise at the output of the receiver. If 
FHSS is present, the jamming power is more efficiently used by 
transmitting all the available power in a limited bandwidth 
which is smaller than the spread spectrum signal bandwidth 
(Ref.4:pp.558]. This strategy is called partial-band inter- 
ference and the fraction of the spread spectrum signal band- 
width which is jammed is denoted by jy. For total jamming 
power J and a FHSS signal bandwidth of W, the barrage noise 
jammer one-sided PSD is N, = J/W over the entire band, while 
the partial-band jammer one-sided PSD is Ni/y = J/yW over a 
bandwidth yW. The total two-sided noise power spectral den- 
sity is y' N,/2+N,/2. The partial-band interference is 
assumed to be present in each branch of MFSK demodulator, and 
the fraction y of the spread spectrum corrupted is assumed to 


be the same for an entire symbol. The fraction y represents 


the probability that partial band interference is present in 
all M branches of the receiver, and l1-y is the probability 
that partial-band interference is not present in all branches 


of the receiver. 


Partial-Band Jamming 


— is 


Total Bandwidth 






Figure 1.2 Barrage and Partial-Band Interference 


C. FADING MULTIPATH CHANNEL. 

In the multipath propagation model for the channel, the 
fading phenomenon is a result of the time variation in the 
Signal phases. Randomly at times, signals can add destruc- 
tively. When that occurs, the resultant received signal can 
be significantly attenuated. At other times, the incoming 
Signals add constructively so that the received signal is 
large. These amplitude variations in the received signal are 
called signal fading and are due to the time variant multipath 
characteristic of the channel[Ref.5:pp.457]. 

Signal fading appears is produced by the signal bouncing 
off mountains and man-made structures as well as atmospheric 
and other interactions such as refraction and diffraction 
processes. The resultant received signal is modeled as a 
complex valued Gaussian process and is assumed to be comprised 
of a direct and many indirect components. When the impulse 
channel response is modeled as a zero mean, complex valued 
Gaussian process, its magnitude envelope at any time is 
Rayleigh distributed. In this case the channel is said to be 
a Rayleigh fading channel. 

In the event that there are signal reflectors in the 
medium, the channel impulse response is no longer be modeled 
as having a zero mean. In this case, its magnitude envelope 
has a Rician distribution and the channel is said to be a 


Rician fading channel. 


The probability density function (pdf) of a Rician faded 


signal is [( Ref.8:pp.108] 





Re az+a? ad 
oh lege Pepe) To) ula) (1) 


where a>0, a* is the power in the direct component, 20’? is the 
power in the diffuse components, I,(a) is the zeroth-order 
modified Bessel function of the first kind, and u(a) is the 
unit step function. For deep fading on the channel, the 
direct communication path is blocked, and as a consequence the 
direct signal component power is zero. In that event, the 
received signal is Rayleigh distributed, and its pdf is 


[Ref.6:pp.108] 


2 
£,(a) = exp [- 5] ula) (2) 


In this thesis, the channel model is assumed to be:1) slow 
fading which implies the signal amplitude remains constant at 
least for the duration of a single hop, and 2) non-frequency 
selective which implies that the hop bandwidth is small 
compared to the coherence bandwidth. This implies that each 
hop frequency can be considered as if it were from an indepen- 
dent channel where the phase and amplitude are relatively 


constant. 


D. BINARY RATIO-STATISTIC COMBING RECEIVER. 

Ratio-statistic combining implies that decision statistic 
for the receiver is a ratio. The general structure of the 
binary ratio-statistic combiner to be considered is depicted 
in Figures 1.3 and 1.4 (Ref.3:pp.26-28]. As illustrated, the 
incoming signal is dehopped and nocoherently demodulated by a 
bandpss filter of bandwidth B=R, (where R, is the hop rate) 
followed by an envelope or square-law detector for each of the 
two BFSK frequencies. The ratio-statistic for each detector 
output of a particular diversity reception is equal to the 
detector output divided by the value of the maximum detector 
eulcpuitee Consequently, one detector output of the ratio- 
statistic receiver for a given diversity reception is always 
unity, and the other detector output for that diversity 


reception is less than unity [Ref.2:pp.146]}. 





Figure 1.3 Diagram of Ratio Statistic Combiner. 
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Figure 1.4 Alternate Implementation of the Ratio-Statistic 


combining receiver. 


The mathematical analysis to develop the probability 
density function for f,(w) was performed by Riley [{Ref.3:pp.27 
-53] where it is shown that Figure 1.3 is equivalent to Figure 
1.4. Furthermore, it is shown that, in Figure 1.4, the random 


variable X, and X, are both normalized to obtain the random 


variable 
x: 
Tk 
where 1= 1 or 2 and q = max(X,,,X,). The random variable 2Z, 


varies between 0 and 1. The random variable Y, varies from -1 


to +1, and the random variable W varies from -L to +L. 


tel 


The probability %ef bit emuer PP Storummamdi versity heandea 


jamming duty factor y is 


L 


P=). ( q)7-y) L-ip. (7) (4) 


Tao 


where P,.(jJ) represents the conditional probability of bit 
error given that j of L hops are jammed. The number of hops 
that are jammed will depend on the duty factor of the jamming 
vy. The total probability of bit error is found by summing 


over all possible values of j. P,.(j) 1s determined by 
0 : 
Pyo(j) = £,(w) dw (5) 


The derivation of the PDF for Y, and W requires conside- 
ration of two cases: when the partial-band interference is 
present and when it is not. The evaluation of f,(w) requires 
an L fold convolution. The noise power at the receiver when 
partial-band interference is present with probability y is 


given by 


o,= (y IN, +N,) B (6) 


eZ 


The noise power at the receiver when no partial-band 


interference is present with probability 1-y is 


of=N_B (7) 


The pdf for the random variable Y, is obtained as follows. 


The dehopped signal is expressed as 


S(t) =/2acos (wt+8) (8) 


where w is the frequency of the dehopped signal, © is the 
phase, (2)°a is the amplitude, and O<t<T, 

If the signal is assumed to be in the upper branch of the 
receiver, then f,,, 1S Rician ditributed and is conditional on 
a being present. The pdf f,, of X,in the lower branch is 
Rayleigh distributed. 

The random variables X, are independent, and the joint 
density function f,4(X,% ,]a) is obtained as the product of 
fy, and fy. Next, appropiate auxiliary variables are 
defined, and the Jacobian of the transformation is obtained. 
Benoiveeeene Conditional pdt f,,(¥,|\a) of Y¥, is obtained by 


integrating over the auxiliary variable. 


a Be 


The unconditional pdf for Y, is 


Evy (Vx) =f Eve (Yqla) £4 (a) da (9) 


where f,(a) represents the pdf of the signal amplitude, which 


is a Rician random variable as given in equation (2). 








The 
diffuse signal-to-noise ratio is defined as 
20° 
Pa 5 (10) 
Ok 
The direct signal to noise ratio is defined as 
a? 
Pr -Z (11) 
Ox 
and the direct to diffuse ratio as 
we 
DDR= (12) 
268 


The PDF for the random variable W as f,(w) is numerically 


evaluated from 


£,,(W) =(fn(¥,) VO ty, (yy (13) 


14 


where the j superscript implies jammed hops and the o super- 
script implies unjammed hops. 
P..(j) 1s obtained by equation (5), and the probability of 


bit error is found by numerically evaluating equation (4). 


1. Unconditional pdf f,,(Y,) for envelope detection 


The conditional pdf for the random variable X,, is 
given by 


(X2,+2a2) 





Ay x 20% V2aX, , 14 
Ex1x (X;|@) = = 2 ee | ; =<) Soe 
O;, 0 
The pdf for the variable X, 1s 
. _ Xak 
2 
Lyon (X24) =—*e ack (15) 
6) 


k 


The probability density function for the random 


variable Y, in the interval of -1 < Y, < 0 is given by 


ibe, 


2 Px 
2 CESS.) 1+€,+(Tey,)* 


rw eee 
vx (Ya) [1+(1+Y¥,)?] (1+8,+(1+Y,) 2] 


(16) 
| Pas lee (1+¥;)* , Px (tt (1+¥,)*) | 


RE ee ee : 1+&,+(P¥Y,) ¢ 


and for 0 < ¥Y,<1 


: p,(1-Y,)? 
PANSY, LS ae 


Ly (Vg) = re er 
Yk K (1+(1-Y,) 2] [I tice) OE ae) =) 


(aa 
a PR 
mabe aye | Sb ea (a si) 





2. Unconditional pdf f,,(Y,) for square law detection. 
The conditional pdf for the random variable xX, is 


given by 


(X2,+2a2) 


¥2a% » , (18) 


2 
Lr (X41) ae met bio 
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The pdf for the random variable xX, is 
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The probability density functions for the random variable Y, 


in the interval 1 < Y¥, < 0 1s given by 


Px 


ab Decoy, 
f,,(Y,) =———--——;+—-- et 
vk (Vy (ZENG) 2a er ae) (20) 
ae (1+y,) (s ‘ Tp Pans) | 
2+ Gree Yee Berg ety 





and for uQe< Y sei 


P,(1-Y;,) 


1 ~ 2=¥yrey (l-Yy) 
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E. FORWARD ERROR CORRECTION CODING. 

Forward error correction (FEC) is a method that employs 
the adding of systematic redundancy and noise averaging at the 
transmit end of a link such that errors caused by the channel 
can be corrected at the receiver by means of a decoding 
algorithm. The binary data source generates information bits 
at R, seconds. These information bits are encoded at a coded 


rate R, and the coded bit energy is 


E_=R E, (22) 
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The encoder output sequence is modulated and transmitted over 
the communication channel. At the receiver, the demodulator 
output is passed through a decoder to recover the original 
binary data. 

FEC can be implemented with either block or convolutional 
codes [Ref.7:pp.417-441]. Convolutional codes are selected 
for this thesis. 

A. J. Viterbi [Ref.8] defines a convolutional encoder as 
a linear finite state machine consisting of a K-stage 
shift register and n linear algebraic function generators. 
The binary input data is shifted along the register b bits at 
a time. 

Each coded symbol carries an average of b/n information 


bits, and the code rate R 1s 


(23) 


ty [o> 


Figure 1.5 depicts the general structure of a convolu- 
tional encoder, where for a given information bit rate the 
required E,/N, for a specific bit error rate with FEC is 
generally less than the E,/N, required without FEC. Coding 
gain is defined as the difference in E,/N, with and without 


FEC |RebS7 pe. 419 | 
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CONSTRAINT LENGTH K=7 





CODED OUTPUT DATA 


Figure 1.5 A Convolutional Encoder. 


The Viterbi decoder is one of the most common decoders for 
convolutional codes because it is the optimum decoding 
algorithm in the sense of maximum likelihood decoding of the 
entire sequence for convolutional codes. The Viterbi decoder 
performance can be approximated through the use of an upper 


bound [Ref.9:pp.407]. For the rate of 1/2, constraint length 


£9 


K=7 convolutional code, the decoded bit error probability is 


bounded by 


Digs < 36D e210 eae el 16 2.8 tee ee ) (24) 


€ 


st 
aa 


where 


D = 2,/P,_(1-P,,. (25) 
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III. METHOD OF RESEARCH 


Computation of the probability of bit error involves a 
numerical evaluation of (12), (33), and (34). Numerical 
evaluations are done for both envelope and square-law 
detection. 

The coded system performance is shown by plotting the 
decoded bit error probability P,. versus the average signal 
bit energy density-to-jamming noise power spectral density 
ratio E,/N;. Several constraints are established in order to 
analyze the performance of the communication system. 

The jamming duty factor y takes the values .001, .05, .25, 
and 1.0 to determine the effects of different amounts of 
partial-band interference on the system. The direct-to- 
diffuse ratio DDR takes the values .01, 10, and 10000 to 
determine the effects of Rayleigh, Rician, and Gaussian 
channels on the system. The bit energy-to-noise ratio E,/N, 
takes the values 13.35 and 16.00 dB to evaluate the effects of 
thermal noise on the system. The bit energy-to-jammer density 
ratio E,/N, takes values from 0 to 40 dB to review the effects 
of jamming against the system and each constraint is evaluated 
for diversity L values from 1 to 4. 

The uncoded and coded performances are compared in order 


to determine the advantages of using FEC with the spread 


zal 


spectrum frequency hopping/noncoherent binary FSK ratio- 
statistic combining receiver in a communication channel 


characterized by fading and partial-band interference. 


a2 


IV. DISCUSSION OF RESULTS 


The channel is modeled for a moderately strong direct-to- 
diffuse signal ratio a’,/20’=10 which is considered Rician, for 
a’ /207=.01 which is considered a Rayleigh faded channel, and 
Gaussian with a’,/20’7=10000. The worst jamming case corres- 
ponds to broadband jamming y=1. 

Performance with FEC in all cases showed an appreciable 
improvement in performance as expected. Comparison of the 
rate 1/2 convolutional code versus uncoded performance shows 
the expected coding gain. The coding gain is based on an 
optimum diversity level, this being the level where the 
probability bit error approches an asymptotic lower limit for 
a specific value of bit energy-to-noise density ratio. 
Uncoded performance showed an asymptotic lower limit of 
P,.=10* for a bit energy-to-noise of 16 dB and P,.=10° for 13.35 
GB; however, decoded performances show an asymptotic lower 
limit of P,.=10° for a bit energy-to-noise density ratio of 16 
db and P,.=10° for 13.35 aB. The worst case occurs at the point 
where the probability of the bit error reaches its maximum. 

Figures 3.1 to 3.4 are the envelope receiver performance 
for partial-band and broadband jamming at a direct-to-diffuse 
ratio of 10 for diversities L=1, 2, 3 and 4. Figure 3.1 shows 


that for bit energy-to-jamming density ratios greater than 10 
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aB, the performance with broadband jamming is better than 25% 
partial-band interference. For bit energy-to-jamming ratios 
greater than 15 dB, again, the performance with broadband 
jamming is better than both 5% and 25% partial-band inter- 
ference. The uncoded system shows performance degradation due 
partial-band interference from 8 to 40 daB, while the coded 
system shows better performance against partial-band inter- 
ference from 15 to 40 dB. For bit energy-to-jamming density 
ratios less than 15 dB, receiver performance improves more for 
partial-band than broadband jamming. 

Figures 3.2, 3.3, and 3.4 show that degradation due to 
partial-band jamming is reduced dramatically for L=2, 3, and 
4. Figure 3.5 shows the worst case for diversities of L=1, 3, 
and 4. It is interesting to note that if we compare the plots 
for diversities of L=1l1 and 4, the better performance 
corresponds to diversity L=4, but between L=3 and 4, diversity 
L=3 performs better than L=4. There are improvements in 
performance when diversity L=3 is used for bit energy-to- 
jamming density ratios above 18 dB. 

Figure 3.6 1s the receiver performance in the presence 
of a Gaussian channel which is degraded significantly for no 
diversity while Figure 3.7 shows how this degradation is 
reduced when diversity is increased. 

Figure 3.8, for a Gaussian channel, we see that the 
performance advantage due to increased diversity is lost for 


bit energy-to-jamming density ratios greater than 26 dB, where 
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the L=1 worst case performance drops below the L=4 worst case 
performance. 

Figures 3.9 and 3.10 show the receiver performance for a 
Rayleigh channel; the worst case performance corresponds to 
broadband jamming and the effects of increasing diversity are 
an important consideration. 

Figure 3.11 shows, for different amounts of channel 
fading, the receiver performance with 5% partial-band inter- 
ference and diversity L=4. Figure 3.12 illustrates the 
effects of broadband jamming. The worst case performance for 
Rayleigh channels corresponds to broadband jamming. 

For E,/N,=13.35 dB, Figures 3.13, 3.14, 3.15, and 3.16 show 
receiver performance for different fractions of partial-band 
and broadband jamming in the presence of a Rician faded 
channel for diversities of L=1, 2, 3 and 4. For no diversity, 
the performance with broadband jamming is better than 5% 
partial-band interference for bit energy-to-jamming density 
ratios from 17 to 30 dB and better than 25% partial-band 
interference for bit energy-to-jamming density ratios lower 
than 30 dB. Partial-band interference degrades the receiver 
performance dramatically versus the receiver performance for 
bit energy-to-noise of 16 dB. When diversity is increased, 
degradation is reduced for diversity L=2 and worst case 
performance occurs for broadband jamming; however, diversity 


greater than L=2 does not improve performance because 
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diversity L=2 has an asymptotic lower limit of P,,=10° while 
for L=4 the asymptotic lower limit is P,.=10%. 

Figure 3.17 shows worst case performance for diversities 
L=1, 3, and 4. The advantage of using diversity is lost for 
low values of bit energy-to-noise density ratios. 

Figure 3.18 shows receiver performance for a Gaussian 
channel, which is better than that for a Rician faded channel 
(GES, Figure" 3486). 

For Gaussian channels, Figure 3.19 shows that performance 
with increased diversity 1s dramatically poorer for all bit 
energy-to-jamming density ratios. 

Figure 3.20 shows both bit energy-to-noise density ratios 
for Rician’ faded’ channels f6r 13.35 and 16 dBag For larger 
values of energy-to-jamming density ratios, receiver perfor- 
mance tends to be controled by thermal noise. 

Receiver performance for L=3 is better than L=1 for a bit 
energy-to-noise density ratio of 16 dB for all bit energy-to- 
jamming density ratios above 17 dB. For a bit energy-to-noise 
density ratio of 13.35 daB, lower orders of diversity perform 
better than higher for all values of bit energy-to-jamming 
density ratios. 

In Figures 3.21, 3.22, and 3.23, the performance of the 
square-law receiver is illustrated for a bit energy-to noise 


ratio of 16 @B and diversities of L=1, 3 and 4 for a Rician 
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channel. For diversities greater than L=1, the worst case 
performance tends to be when the interference is broadband. 
For the envelope detection case no performance improvement is 
obtained for diversities greater than L=3. Figure 3.24 and 
Figure 3.25 show worst case performance for diversities of 
L=1, 3, and 4 for bit energy-to-noise density ratios of 16 dB 
and 13.35 dB, respectively. 

Figures 3.26 and 3.27 show comparisons between both coded 
detections examined for bit energy-to-noise density ratios of 
13.35 and 16 dB. For diversities of L=1 and L=3, the envelope 
detector performs better than the square-law detector for all 
values of bit energy-to-jamming density ratios. For uncoded 
detection, the envelope detector performs better than square- 


law detector for diversities of L=3 or 4. 
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Figure 3.3. Performance of a FFH/NCBFSK/RS Receiver with 


Envelope detection, where E,/N.=16 dB, Direct 


to Diffuse RATIO=10, and L=3. 
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to Diffuse RATIO=10, and the worst case for 


L=1, 3 and 4. 
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Figure 3.7. Performance of a FFH/NCBFSK/RS Receiver with 


Envelope detection, where E,/N.=16 dB, Direct 


to Diffuse RATIO=10000, and L=4. 
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Figure 3.8. Performance of a FFH/NCBFSK/RS Receiver with 


Envelope detection, where E,/N.=16 dB, Direct 


to Diffuse RATIO=10000, and Worst L=1 and 4. 
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Figure 3.9. Performance of a FFH/NCBFSK/RS Receiver with 
Envelope detection, where E,/N.=16 dB, Direct to 


Diffuse RATIO=.01, GAMMA=.05, and L=1,2,3 and 4. 
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Figure 3.10. Performance of a FFH/NCBFSK/RS Receiver with 
Envelope detection, where E,/N.=16 dB, Direct to 


Diffuse RATIO=.01, GAMMA=1, and L=1,2,3 and 4. 
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Figure 3.12. Performance of a FFH/NCBFSK/RS Receiver with 
Envelope detection, where E./N,=16 aB, Different 


Amounts of Fading, GAMMA=1, and L=4. 
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Figure 3.15. Performance of a FFH/NCBFSK/RS Receiver with 


Envelope detection, where E,/N,=13.35 dB, Direct 


to Diffuse RATIO=10, and L=3. 


42 


ENVELOPE DETECTOR Eb/No= 13.35 dB L 


19°@ 


16-1 


18-3 


19-4 


Probability of Bit Error 


19-5 


14-6 
15 


Figure 3.16. 


4 RATIO= 18 





20 25 38 35 486 45 


Eb/NJ 
JAMMING RATIO GAMMA 


*= 001 
o = .05 
-= 25 


-= 1.00 
Performance of a FFH/NCBFSK/RS Receiver with 


Envelope detection, where E,/N.=13.35 dB, Direct 


to Diffuse RATIO=10, and L=4. 


43 


Probability of Bit Error 


ENVELOPE DETECTION Eb/No= 13.35 dB RATIO= 18 


18° 


189-2 


19-3 


19-5 





19-6 
15 
Eb/Nj 
DIVERSITY 
*=L1 o=L3 -=L4 


Figure 3.17. Performance of a FFH/NCBFSK/RS Receiver with 
Envelope detection, where E,/N,=13.35 dB, Direct 
to Diffuse RATIO=10, and the worst case of L=1, 


3 and 4. 
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Figure 3.18. Performance of a FFH/NCBFSK/RS Receiver with 
Envelope detection, where E,/N,=13.35 dB, Direct 
to Diffuse RATIO=10000, L=4 and different amounts 


of interference. 
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Figure 3.19. Performance of a FFH/NCBFSK/RS Receiver with 
Envelope detection, where E,/N.=13.35 dB, Direct 
to Diffuse RATIO=10000, and the worst case L=1 


and 4. 
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Figure 3.20. Comparison of the Performance of a FFH/NCBFSK/RS 


Receiver with Envelope detection, where Direct 
to Diffuse RATIO=10, E,/N.=13.35, and 16 dB and 


worst L=1 and 3. 
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Figure 3.21. Performance of a FFH/NCBFSK/RS Receiver with 


Square Law detection, where E,/N,=16 dB, Direct 


to Diffuse RATIO=10, and L=1. 
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Figure 3.22. Performance of a FFH/NCBFSK/RS Receiver with 
Square Law detection, where E,/N,=16 dB, Direct 


to Diffuse RATIO=10, and L=3. 
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Figure 3.24. Performance of a FFH/NCBFSK/RS Receiver with 
Square Law detection, where E,/N,=16 dB, Direct 
to Diffuse RATIO=10, and the worst case of L=1, 


3, and 4. 
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Figure 3.25. Performance of a FFH/NCBFSK/RS Receiver with 
Square Law detection, where E,/N=13.35 dB, 
Direct to Diffuse RATIO=10, and the worst case 


of L=1, 3, and 4. 
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Figure 3.26. Comparison of the Performance of a FFH/NCBFSK/RS 


Receiver both detection, where E,/N,=16 dB, Direct 


to Diffuse RATIO=10, and the worst case of L=1 


and L=3. 
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Figure 3.27. Comparison of the Performance of a FFH/NCBFSK/RS 
Receiver both detection, where E,/N,=13.35 dB, 
Direct to Diffuse RATIO=10, and the worst case of 


L=1 and L=3. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS. 

This thesis has shown the performance of envelope and 
square-law detection of L-hops/symbol FH/BFSK for worst case 
partial-band jamming. The nonlinear ratio-statistic combining 
technique with forward error correcting coding is an effective 
strategy when combined with diversity to mitigate jamming 
effects. 

Ratio-statistic combining soft decision receivers with 
hOmWatce error mae coding outperforms the uncoded 
receiver for all cases of diversity, noise jamming inter- 
ference, and channel fading. 

The envelope detection optimum diversity level tends to be 
L=3 with a bit energy-to-jamming density ratio of 16 db; this 
corresponds to broadband jamming. 

System performance is dramatically degraded when the 
bit energy-to-noise density ratio is reduced to 13.35 dB. The 
worst case performance for no diversity is better than higher 
orders; however, the worst case performance is better than 
that of the uncoded system under the same conditions. 

The coded performance of square-law detection systems is 


not improved by diversity when jamming is broadband for all 
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channel types. 

Comparision of the performance of envelope detectors and 
square-law detectors show that envelope detection performs 
better than square-law detection for all bit energy-to-jamming 


density ratios considered 


B. RECOMMENDATIONS. 

The FFHSS/noncoherent BSFK receiver with ratio-statistic 
combining and forward error correcting coding has been 
evaluated for partial-band jamming and for different types of 
faded channels. This analysis should be extended to the 
general M-ary case. In addition, the ratio-statistic system 
should be simulated using some of the available communication 


Simulation software packages. 


JO 


APPENDIX 


MATLAB SOURCE CODES 


% The following source codes were developed by 

% JOHN RILEY, 1990 and this codes were modified 
% to implement the 1/2 VITERBI decoder algorithm. 
% ENVMAIN.M Envelope Detection. 

% This is the main source code that generates 

% the decoded Pbhe for the system with Ratio 

% Statistic combiner and convolutional VITERBI 

% decoder of rate 1/2. 

$ 

me 1/2; 

Real1O >>) INPUT TEHE VALUE OF DIRECT TO DIFFUSE RATIO ’; 
n= 150; 

Eb = 10*(xX.xXxX);’ INPUT THE VALUE OF BIT ENERGY’ 

Ec = r*Eb 

No = 1.0; 


GAMMA =’ INPUT THE VALUE OF FRACTION JAMMING ” ; 
i=" SNPUT THE DESERE SYMBOL DIVERSITY ‘; 
Eh = Ec/L; 
DIFFUSE = Eh/ (RATIO+1); 
DIRECT = DIFFUSE*RATIO; 
RHOn = DIRECT; 
Miie= DLP EUSE > 
ROr XK S121 
clear Pb 
EDN} = (XxX~-1)*2; 
Nj} = 10*(-EbNjJ/10) *Ec; 
Nt = 1 + Nj/GAMMA; 
RHO = DIRECT/Nt; 
hi = DEFFUSE/NE > 
envpbe; 
PPbe(xx) = Pb; 
Meee —s10* jlogi0(Eb/N}) ; 
end 
save PPbe; 
save X; 
clear 
semllogy (X,PPbe) ; 
title(’Eb/No= dB GAMMA= L= RATIO= 
Xlabel(’Eb/N}j’); 
Vilanema Peobabi lity of Bit Error’); 


a | 


% ENVPBE.M for rate code of 1/2 
% This program will determine the probability of error 
% as in equation 5 of Keller’s paper. It is the 
% algorithm for the 1/2 VITERBI decoder for envelope 
% detection. 
Pbhee = 0; 
for j = 02 
convfad; 
simperr; 
Pbc (j+1)=Pbej* (fact(L) /(fact(j)*fact(L-j) ) ) *GAMMA*74* (1-GAMMA 
) (Eaape 
end 


Pbee=sum(Pbc) 
d = 2*sqrt((1-Pbee) *Pbee) ; 
valor (xx) =d 


Gq = —36*~(alaero; 
Ww =e 21 dane 
ee = M4eaed)-14- 
Af = 11633*(da)*16; 


Pb= .5*(qq+wwtee+Af) ; 
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FADNOINT.M Envelope detection 

This program will sample the probability density 
functionm of the variable y, it generates the NON 
JAMMED HOPS sequence and put it into and array f(y). 


HP WP oP oP ot? 


lear f 


II 
mMAMOANNFO 


ae 


> 


erie? x 
oo 

-1+2/n; 

2 ereyoe cl iaty) 2) (itl ity) ~2+xXin) ); 
Roem, i+ (1l+yV) 2+XiEn) ; 

(ive / (Ln ty ) 72) 4 
RHGma(iae ity) 2) 7 (1 +XInt(1+y) *2) ; 

= C*exp(-D) *(1+E* (XIn+kK) ) ; 


lou i tt te il 


end 
Gor xX = olagall 

Ke maeilis 

Zeya 

2 eye) Met Loy) 2) s (lt (ly) o2*(1+Xin) )) ; 
RHGwet yy) 2) lat 1 -y) 2* (1+xXin)); 

tele GX it) eee yea 2) 

RHOn# (1+ (1-y) *2) / (1+ (XIn+1) *(1l-y) *2) ; 
(xtn) = C*exp(-D) *(1+E* (XIn+K) ); 


II 


loi Ww we ue 


1 
Z 
y, 
E 
D 
E 
K 
f 


end 
Ve— =te2/length(£):1—2/length(f£) ; 
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end 
Lore 


end 


il 
mMAMOAN NFO 


ra 


il 
mMAMOAN NP 


-~ 


FADWINT.N Envelope detection 

This program will sample the probability density 
function of the variable y, it generates the JAMMED 
HOPS sequence and put it into and array f(y). 


lear f 


> eae | 
=lt27 i, 
2*( ity) /((1+(14+y) 22) 2G Gy 2 ae ey, 
RHOi/ (1+(1+y) *2+XIi); 
(1+y) *2/(1+XI1i+(1+y) *2); 
RHOi* (1+(1+y) *2) / (1+XIi+ (lity) *2) ; 
= C*exp(-D) *(1+E*(XIi+K) ); 


lulu it wt ws 


s<n+1 

xX le 

z/n; 

2* (1=y)/ (1 (1 =y) ~2)) * ae een ea); 
RHO1* (1-y) *2/ (1+(1-y) *2* (1+XI1) ) ; 

1/ (iGo it) elas) 2 

RHOi* (1+(1-y) *2) /(1+(XIit1) *(1-y) *2); 
+n) = C*exp(-D) *(1+E* (XIi+kK) ) ; 


“mi lt uw te ue il 


Vv = -1:2/length(f):1-2/length(f) ; 
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2 / 
j== 


if 
elseif 
elseif 
elseif 
elseif 
elseif 
elseif 
else 

end 

fadnol 
if (L- 
elseif 
elseif 
elseif 
elseif 
elseif 
elseif 


else 


end 
fw 


CONVFAD.M ENVELOPE DETECTION 

This rutine generates the convolution of the Folded 
convolutions for the jammed and non jammed HOPS. 

fw represents the numerical evaluated PROBABILITY 
DENSITY FUNCTION for the random variable W. 


e 
(length(f)-1); 
fwi = (1 ]; 
j==1 
fwl = f; 
j== 
cw "= Com (tL, f); 
j== 
Ewim—meom (COnvV(£,f),£) ; 
== 
iWiwo—comu( CONV (convy(£,f),f£),£); 
== 
Eve Cony (Comv( Conv (conv(£,f£f),£),£),£); 
j== 
Evin SCO conv Conv(conv(conv(f,f),f£),£),£),£); 
fwi = error; 
nt 
j) == 
fwn = {( 1 ]; 
Cos] ea 
fwn = f; 
a2 
mum = cConvi tf, £) > 
CS) se 
Pvneaeecony(cCOonv(L,£) ,£) ; 
(L-j)==4 
fwn = conv(conv(conv(f,f),f£),£); 
ea) a= 
Ewe cony (Cony conv(conv(f,£),f£),f£),£); 
(3) 6 
[Vie wean | CONV (COnyY(COonvy(Conv(£,f£),f£),£),£),£); 
fwn = error; 


h* (L-1) *conv(fwi, fwn) ; 


aj 


HO HAP AP DP OW AP AO 


SIMPERR.M 

This program uses Simpson’s rule to find the area 

under the array fw for y less than zero. This will 
be the probability of bit error for a given # HOPS 
jammed, Pbhe(j) in a receiver using Ratio-Statistic 
Combining. The result is in Pbej. 


areal=0; 
Pbej = 0; 
2: (lengen(i~.)-—1))7 2-1 
if round(s/2) == s/2 

areal = areal + 4*fw(s); 
else 

areal = areal + 2*fw(s); 
end 


areal + fw(1) +fw((length(fw)-1)/2); 
areal*h/3; 
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SQMAIN.M SQUARE LAW DETECTION 
This is the main source code that generates 
the decoded Pbhe for the system with Ratio 
Statistic combiner and convolutional VITERBI 
decoder of rate 1/2. 


FY AP HW AP AP AP OP 


= 1/2; 
RATIO = ‘INPUT THE VALUE OF DIRECT TO DIFFUSE RATIO ’; 
n= 150; 
Eb = 10*(x.xXxX);’ INPUT THE VALUE OF BIT ENERGY’ 
Ec = r*Eb 
No = 1.0; 


GAMMA =’ INPUT THE VALUE OF FRACTION JAMMING “7 ; 
fio INPUT STHE DESTRE SYMBOL DIVERSITY ‘’; 
Eh = Ec/L; 
DIFFUSE = Eh/(RATIO#1) ; 
DIRECT = DIFFUSE*RATIO; 
RHOn = DIRECT; 
XIn = DIFFUSE; 
for xx =1:21 
clear Pb 
EbDNJ = (xXx-1) *2; 
Nj = 10*(-EbDNjJ/10) *Ec; 
Nt = 1 + NjJ/GAMMA; 
RHOi = DIRECT/Nt; 
XIi = DIFFUSE/Nt; 
sqpbe; 
PPbe(xx) = Pb; 
X(xx) = 10 * logl0(Eb/N3); 
end 
save PPbe; 
save X; 
clear 
semilogy(X,PPbe) ; 
title(’Eb/No= dB GAMMA= L= RATIO= 
walabel(“Bb/ Na”); 
Vaabel~ Eeobabtlity of Bit Error’); 
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SQCORRE.M for rate code of 1/2 

This program will determine the probability of error 
as in equation 5 of Keller’s paper. It is the 
algorithm for the 1/2 VITERBI DECODER for the square 
law detection. 


OP AP AP OP OP CNP 


Phee = 0; 
for 3 *="02% 
sqcovfad; 
sSimperr; 
Pbc (j+1)=Pbej* (fact (L) /(fact(j) *fact(L-j) ) ) *GAMMA*4*(1-GAMMA 
ie Gi=2))c, 
end 


Pbee=sum(Pbc) 
d = 2*sqrt((1-Pbee) *Pbee) ; 
valor (xx) =d 


qq = 36*(d)%*10; 
Cs 2s sl) OP ¢ 
Ce = lear (ciianiean 
Af = ll6sexrd)oter 


Pb= .5*(qqt+twwteetAf) ; 
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dO AP AO DO OP 


more 


end 
For 


end 


Vv = 


SQFADNOI.M Square law detection 

This program will sample the probability density 
function of the variable y ,it generates NON JAMMED 
HOPS sequence and put it into and array f(y). 


ou 1 

-1+2/n; 

ee ee2 ty) *(xint2ty))); 
(1+y) /(XIn+2+y) ; 

=Rinom, (,in+2+y) ; 
RHOn+~(2+y) / (Xint2+y) ; 

= Cxexp(E) *(1+D* (XIn+K) ) ; 


sot tt t tf uo 


al 


sn+1 

a ne 

Z/ ine 

a ee ene (1-V)) * (2-y)); 
RHOn* (1-y) / (2-y+XIn*(1-y)) ; 
Deny eee iy) js 

RHOnAK sy) / (2mytXIn* (1-y)y); 
(xtn) = C*exp(-D}*(1+E* (XIn+kK) ) ; 


=li2y7 benden nj: £—2 / length (f£).; 
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SQFADWIN.M Square law detection 

This program will sample the probability density 
function of the variable y ,it generates the JAMMED 
HOPS sequence and put it into and array f(y). 


oO ot? 


lear f 


ll 
mA MOON NFA 


~7“-~ 


hor 1X 
x - 1; 

-1+2z/n; 

(1/ ( (2+y) *(XI1+2+y))); 
(1+y) 7 Cyt tore oe, 

-RHOi/ (XIi+2+y) ; 
RHOi*(2+y) /(XIi+2+y) ; 

= C*exp(E) *(1+D*(XIi+K) ); 


x“ ui ti wt te ues 


ball 


end 
for x 


mAMOAN NF 


~—, 


snl 

x = a 

VAIO 

1/ (( 25y Pat aS i= (2 Se, 
RHOa* (1 sy) (2 ava (lave, 
l/ (2sy +xXit (rae 

RHO1 * (Za (2 a ee, 
+n) = C*exp(-D)*(1+E*(XIi+k) ); 


x“ i uw wt wt ot 


end 
V = -1:2/length(f) :1-2/length(f) ; 
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% SQCOVFAD .M SQUARE LAW DETECTION 
% The rutine generates the convolution of the Folded 
% convolutions for jammed and nonjammed HOPS. 
% fw represents the numerical evaluated PROBABILITY 
% DENSITY FUNCTION for the random variable W. 
% 
sqfadwin 
ho 2/1 Lengen (rt) -l), 
if j== 
fwi = { 1]; 
elseif j== 
fwl = f; 
elseif j== 
fwi = conv(f,f); 
elseif j== 
Gwe mcomy Ceonv (f,£) ,£) ; 
elseif j== 
twimowocony (aony (conv(f£,f£),f£),£) ; 
elseif j==5 
tVvin—  Conmuoeny (conv (Conv(f,f),f),£),f£) > 
elseif j== 
Ryo HConvwe coma coh. conv (conv(f,f),£),£),f),f); 
else 
fwl = error; 
end 
sqfadnol 
if (L-j)==0 
Ewie— { 1) 
elseif (L-j)==1 
fwn = f; 
elseif (L-j)==2 
fwn = conv(f,f); 
elseif (L-j)==3 
nw — GConvy(conv(f,f), f). 
elseif (L-j)==4 
EWiemeaconvy (convy(conv(f,f),f),f); 
elseif (L-})==5 
fwn = conv(conv(conv(conv(f,f),f£),£),£); 
elseif (L-j) ==6 
fii Com (comy (conv ( CONV (COnV(E,£),£),£),£),£); 
else 
fwn = error; 
end 


Eve ene(b-l reconv (fwi,fwn) ; 
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